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TRIBES AND CLADES WITHIN APIACEAE 
SUBFAMILY APIOIDEAE: THE CONTRIBUTION OF 
MOLECULAR DATA 


S. R. DowwiE*T, G. M. PLUNKETT}, M. F. Watsons, K. SPALIK**, 
D. S. KArz-DowwniE*, C. M. VaLIEJO-ROMANTT, E. I. TERENTIEVATT, 
A. V. TnorrskvTT, B.-Y. LEE*, J. LAHHAMS$8 & A. EL-OQLAHŞŞ 


Phylogenetic analyses of chloroplast gene (rbcL, matK ), intron (rp/16, rps16, rpoC1) 
and nuclear ribosomal DNA internal transcribed spacer (ITS) sequences and 
chloroplast DNA restriction sites, with supplementary data from variation in size of 
the chloroplast genome inverted repeat, have been used to elucidate major clades 
within Apiaceae (Umbelliferae) subfamily Apioideae Drude. This paper summarizes the 
results of previously published molecular cladistic analyses and presents a provisional 
classification of the subfamily based on taxonomic congruence among the data sets. 
Ten tribes (Aciphylleae M. F. Watson & S. R. Downie, Bupleureae Spreng., Careae 
Baill., Echinophoreae Benth., Heteromorpheae M. F. Watson & S. R. Downie, 
Oenantheae Dumort., Pleurospermeae M. F. Watson & S. R. Downie, Pyramidoptereae 
Boiss., Scandiceae Spreng. and Smyrnieae Spreng.) are erected or confirmed as 
monophyletic, with Scandiceae comprising subtribes Daucinae Dumort., Scandicinae 
Tausch and Torilidinae Dumort. Seven additional clades are also recognized but have 
yet to be treated formally, and at least 23 genera examined to date are of dubious 
tribal or clade placement. The utility of these different molecular markers for 
phylogenetic inference in Apioideae is compared based on maximum parsimony 
analyses of subsets of previously published molecular data sets. Of the six loci 
sequenced, the ITS region is seen to be evolving most rapidly and rbcL is the most 
conservative. Intermediate in rate of evolution are matK and the three chloroplast 
introns; with rp/16 and rps16 evolving slightly faster than matK or rpoC1. The analysis 
of restriction sites, however, provided 2-4 times more parsimony informative 
characters than any single DNA locus sequenced, with estimates of divergence just 
slightly lower than that of the ITS region. The trees obtained from separate analyses of 
these reduced data sets are consistent with regard to the major clades inferred and the 
relationships among them. Similar phylogenies are obtained by combining data or 
combining trees, representing the supermatrix and supertree approaches to 
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phylogenetic analysis, respectively. The inferred relationship among the tribes and 
informally recognized major clades within Apioideae is presented. 


Keywords. Chloroplast intron, internal transcribed spacer, molecular systematics, 
nuclear rDNA, phylogeny, suprageneric classification, Umbelliferae. 


INTRODUCTION 


‘Although Umbelliferae were the first family of flowering plants to achieve general 
recognition, after nearly three and a quarter centuries of successive and multinational 
effort, considerable disagreement still exists as to the proper delimitation of the family 
and even more uncertainty prevails as to its natural subdivisions and the criteria on 
which they should be erected. Clearly continued acquisition of new information and 
re-examination, refinement, and re-evaluation of accumulated evidence are urgently to 
be welcomed.’ 

L. Constance (1971) 


The higher-level relationships within Apiaceae (Umbelliferae) have been traditionally 
difficult to resolve, particularly within its largest subfamily Apioideae Drude. 
Comparison of the accounts of, for example, Koch (1824), de Candolle (1830), 
Bentham (1867), Boissier (1872), Drude (1898), Calestani (1905) and Koso- 
Poljansky (1916), which were erected largely on the basis of fruit morphology and 
anatomy, shows widely diverging opinions on the definition and composition of its 
tribes and subtribes. While such focus on fruit structure has been rejected by many 
(Heywood, 1971b, 1978a; Theobald, 1971; Davis, 1972; Cronquist, 1982; Hedge 
et al., 1987; Shneyer et al., 1992, 1995), the highly criticized century-old system of 
Drude (1898), or some modification thereof (Pimenov & Leonov, 1993), remains 
the most commonly used treatment. As such, many of the authors who published 
papers as a result of two international symposia on the family (Heywood, 1971a; 
Cauwet-Marc & Carbonnier, 1982) presented their findings in the framework of 
Drude’s system. 

Drude (1898) divided Apiaceae into three subfamilies (Apioideae, Hydrocotyloideae 
Link and Saniculoideae Burnett), recognizing 8 tribes and 10 subtribes within 
Apioideae. Molecular phylogenetic studies, however, while confirming the monophyly 
of subfamily Apioideae, have shown that many of its tribes and subtribes are not 
monophyletic (Downie & Katz-Downie, 1996; Kondo et al., 1996; Plunkett et al., 
1996a,b, 1997; Downie et al., 1996, 1998, 2000a,b,c; Valiejo-Roman et al., 1998; 
Katz-Downie et al., 1999; Plunkett & Downie, 1999). The apioid umbellifers display 
a remarkable array of morphological and anatomical modifications of their fruits, 
some of which are likely adaptations for various modes of seed dispersal (Jury, 
1986). Not surprisingly, these characters are prone to convergence and their almost 
exclusive use to delimit suprageneric taxa has confounded attempts to identify mono- 
phyletic groups. Our goal over the past six years has been to resolve the higher-level 
relationships within subfamily Apioideae. This is necessary in order to provide 
the framework for lower-level revisions, as well as to interpret patterns in the evolu- 
tion of the many phytochemical, anatomical, cytological, morphological and 
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palynological characters available for the group (such as those summarized in 
Heywood, 1971a, and Cauwet-Marc & Carbonnier, 1982). Although these characters 
have been surveyed widely, trends in their evolution and their reliability in demarcat- 
ing taxonomic groups have rarely been considered outside of the framework of 
Drude's system (for exceptions see Plunkett et al., 1996b, and Katz-Downie et al., 
1999). Eventually, the synthesis of these newly acquired molecular data with a 
re-evaluation of existing evidence will culminate in the production of a modern 
classification for the subfamily. 

This paper has three major objectives. First, we summarize the results of a variety 
of previously published macromolecular analyses as they apply to clarifying the 
evolutionary history of subfamily Apioideae. To date, major clades within the 
subfamily have been elucidated on the basis of phylogenetic analyses of chloroplast 
gene (rbcL, matK ), intron (rp/16, rps16, rpoC1) and nuclear ribosomal DNA ITS 
sequences and chloroplast DNA (cpDNA ) restriction sites, with supplementary data 
obtained from variation in size and position of the chloroplast genome inverted 
repeat. Based on taxonomic congruence among these previously published data sets, 
a provisional classification of the subfamily is proposed. Second, by focusing on 
subsets of these molecular data, we compare the utility of different molecular markers 
for phylogenetic inference within the subfamily. Third, based on separate and com- 
bined analyses of these subset data (or the trees derived from them) and the results 
of our earlier publications using non-reduced data sets, the phylogenetic relationships 
among the major clades comprising Apioideae are inferred. While our studies are 
still very much in progress, the results obtained to date provide the necessary frame- 
work and explicit phylogenetic hypotheses from which future revisions and evolution- 
ary studies can proceed. Authorities for names at rank genus and below are given 
in Tables 2-5; authorities for suprageneric taxa and others not in these tables are 
given when first used. 


SuMMARY OF PUBLISHED MACROMOLECULAR ANALYSES 


The earliest comparative studies of nucleic acids in Apiaceae were the DNA-DNA 
hybridization experiments of Valiejo-Roman et al. (1979, 1982) and Antonov et al. 
(1988; Table 1). The groups obtained on the basis of these hybridization data, while 
supporting the distinctiveness of subfamily Apioideae and its separation from subfam- 
ily Saniculoideae, stood in stark contrast to any existing system of classification 
available for Apioideae. DNA distances among the resolved apioid groups, however, 
were sufficiently low and the sampling too sparse to propose an alternative tribal 
treatment. Relationships within the family were also inferred using immunological 
comparisons of seed storage proteins. The earliest of these studies (Pickering & 
Fairbrothers, 1970) supported Drude's division of the family into three subfamilies; 
it also indicated that Apioideae are more similar serologically to Saniculoideae than 
to Hydrocotyloideae. Further analyses by Pickering & Fairbrothers (1971) confirmed 
five serological groupings for nine genera within Apioideae, corresponding to Drude's 


TABLE 1. A summary of systematic studies of Apiaceae subfamily Apioideae incorporating macromolecular data 


DNA-DNA Systematic RFLPs or restriction rpll6  rpsl6  rpoCl Chloroplast IR. 
Reference hybridization serology site mapping ITS intron intron intron matK rbcL structure 


Antonov et al. (1988) e 

Debonte et al. (1984) e 

Downie & Katz-Downie (1996) e 

Downie & Katz-Downie (1999) e 

Downie et al. (1996) e 

Downie et al. (1998) e e 

Downie et al. (20002) e 

Downie et al. (2000b) e e 

Downie et al. (2000c) e e 

Katz-Downie et al. (1999) e 

Kondo et al. (1996) e 
Lee & Downie (1999) e 

Lee & Downie (2000) e e 

Lee & Rasmussen (1998) e 

Mitchell et al. (1998) e 

Pickering & Fairbrothers (1970) ° 

Pickering & Fairbrothers (1971) 
Pimenov et al. (1999) e e 

Plunkett & Downie (1999) e 

Plunkett & Downie (2000) e 
Plunkett et al. (19962) e 

Plunkett et al. (1996b) ° 

Plunkett et al. (1997) e e 

Radford et al. (2001) e 

Shneyer et al. (1991) e 

Shneyer et al. (1992) e 


poe 
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TABLE 1. (Continued.) 


Reference 


Shneyer et al. (1995) 

Soltis & Kuzoff (1993) 
Soltis & Novak (1997) 
Valiejo-Roman et al. (1979) 
Valiejo-Roman et al. (1982) 
Valiejo-Roman et al. (1998) 
Vivek & Simon (1998) 
Vivek & Simon (1999) 
Watanabe et al. (1998) 


DNA-DNA Systematic RFLPs or restriction 


hybridization serology site mapping ITS 
e 
e 
e 
e 
e 
e 
e 
e 
e 


rpll6 
intron 


rps16 
intron 


rpoCl 
intron matK rbcL 


Chloroplast IR 
structure 
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tribes Dauceae W. D. J. Koch, Scandiceae Spreng., Coriandreae W. D. J. Koch, 
Apieae Drude and Peucedaneae Dumort. Representatives of the last three tribes had 
the greatest protein similarity, whereas those from Dauceae and Scandiceae were 
each serologically distinct. These serological data also supported Drude's division 
of Peucedaneae into three subtribes: Angelicinae Tausch, Peucedaninae Tausch and 
Tordyliinae Drude. In contrast, the analyses of Shneyer et al. (1991, 1992, 1995), 
also based on immunochemical reactions of seed storage proteins but with greater 
sampling in tribes Peucedaneae and Smyrnieae Spreng., clearly indicated that many 
of the tribes and subtribes recognized by Drude do not form monophyletic groups. 
Subtribes Angelicinae and Peucedaninae were determined to be serologically not very 
well separated, whereas subtribe Tordyliinae was maintained as distinct. 

The past several years has seen a flurry of molecular systematic activity on subfam- 
ily Apioideae ( Table 1). Many of these studies considered DNA sequences from both 
coding and non-coding loci of nuclear and chloroplast genomes, with fewer studies 
incorporating cpDNA RFLP or structural rearrangement data. Considering only 
GenBank accessions of nucleotide sequence data, and excluding Apium graveolens, 
Daucus carota and Petroselinum crispum (because much of the data available for 
these economically important apioid species are the result of non-systematic studies), 
there has been a fivefold increase in the number of submitted sequences in over just 
the past two and a half years (from c.300 in March 1998 to c.1500 in September 
2000), primarily the result of these comparative phylogenetic analyses. Considering 
all molecular phylogenetic studies cited in Table 1, a total of 195 apioid genera and 
450 species have been examined to date, representing almost half of the 404 genera 
within the subfamily recognized by Pimenov & Leonov (1993). 


Subfamilial relationships 


The results of earlier molecular systematic investigations have confirmed the mono- 
phyly of subfamily Apioideae and revealed its sister group to be subfamily 
Saniculoideae (Downie & Katz-Downie, 1996, 1999; Plunkett et al, 1996a, 1997; 
Downie et al., 1998, 2000b). Saniculoideae are also monophyletic upon the removal 
of Lagoecia (Plunkett et al., 1996b; Downie et al., 2000c). However, the circumscrip- 
tion of the largely herbaceous subfamily Saniculoideae is confounded by its affinity 
to the woody apioid genera Steganotaenia Hochst. and Polemanniopsis B. L. Burtt 
(Fig. 1; Downie & Katz-Downie, 1999). The incorporation of Steganotaenia and 
Polemanniopsis into an expanded Saniculoideae (and their treatment as a separate 
tribe or subtribe) or their recognition as a new subfamily adjacent to Saniculoideae 
is subject to further investigation. Hydrocotyloideae are polyphyletic, with some 
genera allied with Apioideae plus Saniculoideae and other genera with Araliaceae 
(Plunkett et al., 1996a, 1997; Downie et al., 1998, 2000a). Upon the exclusion of the 
latter group (i.e. the ‘araliaceous hydrocotyloids’ sensu Plunkett et al., 1997; Fig. 1), 
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Apioideae 


o 
Saniculoideae S 
S 
a 
< 
Steganotaenia & £ 
Polemanniopsis oO 


Azorella clade 

(Azorella, Bolax & Eremocharis 
and possibly Bowlesia & 
Klotzschia) 


Araliaceous 


Hydrocotyloids 

(Centella, Didiscus, Hydrocotyle, 
Micropleura, Spananthe, & Xanthosia) 
& Araliaceae 


FIG. 1. Summary of relationships within core Apiaceae as revealed by phylogenetic analyses 
of molecular data (Plunkett et al., 1996a, 1997; Downie & Katz-Downie, 1999; Downie et al., 
1998, 2000b). Members of the Azorella clade and the araliaceous hydrocotyloids (Didiscus 
DC. ex Hook. = Trachymene Rudge) are treated by Pimenov & Leonov (1993) in Apiaceae 
subfamily Hydrocotyloideae; the araliaceous hydrocotyloids do not comprise a monophyletic 
group (Plunkett et al, 1997; Downie & Katz-Downie, 1999; Downie et al., 2000b). 
Steganotaenia and Polemanniopsis (Apioideae and incertae sedis, respectively; Pimenov & 
Leonov, 1993), form a clade sister to Saniculoideae and, pending further study, may belong 
within this subfamily (Downie & Katz-Downie, 1999). 


Apiaceae comprise a clade. Sister to Apioideae plus Saniculoideae, and collectively 
forming the ‘core Apiaceae (Plunkett et aL, 1997)’, are three genera previously 
attributable to Hydrocotyloideae (Azorella Lam., Bolax Comm. ex Juss. and 
Eremocharis Phil.). We have named this group the Azorella clade (Fig. 1; Downie 
et al, 2000b), fully aware that the subfamilial name Azorelloideae, erected by 
Cerceau-Larrival (1962), is invalid. Bowlesia Ruiz. & Pav. and Klotzschia Cham. are 
closely related to the Azorella clade (Plunkett et al., 1996a; Downie et al., 1998; 
Downie & Katz-Downie, 1999); all five genera, however, have yet to be analysed 
simultaneously. 


Tribes of Apioideae 


Within Apiaceae subfamily Apioideae, ten tribes and three subtribes have been erected 
or supported as monophyletic as a result of previously published phylogenetic analy- 
ses of molecular data (Table 2). These include tribes Aciphylleae M. F. Watson 
& S. R. Downie, Bupleureae Spreng., Careae Baill., Echinophoreae Benth., 
Heteromorpheae M. F. Watson & S. R. Downie, Oenantheae Dumort., 
Pleurospermeae M. F. Watson & S. R. Downie, Pyramidoptereae Boiss., Scandiceae 


TABLE 2. Generic composition of tribes and subtribes within Apiaceae subfamily Apioideae supported or erected on the basis of phylogenetic 
analyses of molecular data. Previously designated names for these clades (according to Downie et al., 1998, 2000a,b,c, and Katz-Downie et al., 
1999) are provided in footnotes. Genera assigned to more than one suprageneric category may not necessarily be monophyletic within that group 


(see original papers for discussion) 


Aciphylleae M. F. Watson & 
S. R. Downie* 


Heteromorpheae M. F. Watson & | Pleurospermeae M. F. Watson & 
S. R. Downie* S. R. Downie? 


Aciphylla J. R. Forst. & 
G. Forst. 

Anisotome Hook. f. 

Gingidia J. W. Dawson 

Lignocarpa J. W. Dawson 

Scandia J. W. Dawson 


Bupleureae Spreng.” 
Bupleurum L. 


Careae Baill.* 
Aegokeras Raf. 
Aegopodium L. 

Carum L. 
Cyclospermum Lag. 
Falcaria Fabr. 
Fuernrohria K. Koch 
Grammosciadium DC. 
Rhabdosciadium Boiss. 


Echinophoreae Benth.‘ 
Dicyclophora Boiss. 
Echinophora L. 
Pycnocycla Lindl. 


Anginon Raf. 

Dracosciadium Hilliard & 
B. L. Burtt 

Glia Sond. 

Heteromorpha Cham. & 
Schltdl. 

Polemannia Eckl. & Zeyh. 


Oenantheae Dumort.f 
Berula W. D. J. Koch 
Cicuta L. 

Cryptotaenia DC. 

Cynosciadium DC. 

Helosciadium W. D. J. Koch 

Lilaeopsis Greene 

Limnosciadium Mathias & 
Constance 

Neogoezia Hemsl. 

Oenanthe L. 

Oxypolis Raf. 

Perideridia Rchb. 

Ptilimnium Raf. 

Sium L. 


Aulacospermum Ledeb. 
Eleutherospermum K. Koch 
Molopospermum 

W. D. J. Koch 
Physospermum Cuss. 
Pleurospermum Hoffm. 


Pyramidoptereae Boiss." 


Bunium L. 

Crithmum L. 

Elaeosticta Fenzl 

Lagoecia L. 

Oedibasis Koso-Pol. 
Pyramidoptera Boiss. 
Scaligeria DC. 
Trachyspermum Link (in part) 


Scandiceae Spreng) 


Daucinae Dumort. 
Agrocharis Hochst. 
Ammodaucus Coss. 
& Durieu 
Cuminum L. 


Daucus L. 

Laser Borkh. ex P. Gaertn., 
B. Mey. & Schreb. 

Laserpitium L. 

Melanoselinum Hoffm. 

Monizia Lowe 

Orlaya Hoffm. 

Pachyctenium Maire & 
Pampan. 

Polylophium Boiss. 

Pseudorlaya (Murb.) Murb. 

Thapsia L. 


Scandicinae Tausch 


Anthriscus Pers. 

Athamanta L. (in part) 
Balansaea Boiss. & Reut. 
Chaerophyllum L. 
Conopodium W. D. J. Koch 
Geocaryum Coss. 

Kozlovia Lipsky 

Krasnovia Schischk. 
Myrrhis Mill. 

Myrrhoides Fabr. 
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TABLE 2. (Continued.) 


Neoconopodium ( Koso-Pol.) 
Pimenov & Kljuykov 

Osmorhiza Raf. 

Scandix L. 

Sphallerocarpus DC. 


è Aciphylla clade (in part). 
> Bupleurum clade. 
* Aegopodium clade. 


3 Nirarathamnos Balf. f. and Rughidia ined., find affinities with this group. 


* Heteromorpha clade. 
* Oenanthe clade. 


Tinguarra Parl. 
Todaroa Parl. 
Torilidinae Dumort. 
Astrodaucus Drude 
Caucalis L. 


Chaetosciadium Boiss. 
Glochidotheca Fenzl 
Lisaea Boiss. 


Szovitsia Fisch. & C. A. Mey. 


Torilis Adans. 


* Physospermum clade. Molopospermum is tentatively included here based on Shneyer et al. (1992). 


^ Crithmum clade. 


Turgenia Hoffm. 
Yabea Koso-Pol. 


Smyrnieae Spreng.* 
Lecokia DC. 
Smyrnium L. 


i Daucus clade. Artedia L. falls in either Daucinae or Torilidinae, depending upon the study (Lee & Downie, 2000). Ferula (in part; Downie 
et al., 2000c) and Glaucosciadium B.L. Burtt & P.H. Davis (K. Spalik, unpubl. data) may each constitute additional separate lineages within 


the tribe. 
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TABLE 3. Major clades within Apiaceae subfamily Apioideae that have yet to be treated formally. Genera assigned to a particular clade may 
not necessarily be monophyletic within that clade (see original papers for discussion) 


Angelica Clade* 
Aethusa L. 
Aletes J. M. Coult. & Rose 
Angelica L. 
Carlesia Dunn. 
Chamaele Miq. 
Chymsydia Albov 
Cnidiocarpa Pimenov 
Cnidium Cuss. (in part) 
Coelopleurum Ledeb. 
Cortia DC. 
Cymopterus Raf. 
Dystaenia Kitag. 
Endressia J. Gay 
Exoacantha Labill. 
Glehnia F. Schmidt ex Miq. 
Grafia Rchb. 
Harbouria J. M. Coult. & Rose 
Imperatoria L. 
Karatavia Pimenov & Lavrova 
Libanotis Haller ex Zinn 
Ligusticum L. (in part) 
Lomatium Raf. 
Meum Hill 
Musineon Raf. 
Neoparrya Mathias 
Notopterygium H. Boissieu 
Oreonana Jeps. 
Oreoxis Raf. 
Orogenia S. Watson 


Paraligusticum V. N. Tikhom. 

Peucedanum L. 

Phlojodicarpus Turcz. ex Ledeb. 

Podistera S. Watson 

Polytaenia DC. 

Pseudocymopterus J. M. Coult. & Rose 

Pteryxia (Nutt. ex Torr. & A. Gray) 
J. M. Coult. & Rose 

Selinum L. 

Seseli L. 

Shoshonea Evert & Constance 

Spermolepis Raf. 

Sphenosciadium A. Gray 

Taenidia (Torr. & A. Gray) Drude 

Tauschia Schltdl. 

Thaspium Nutt. 

Tommasinia Bertol. 

Xanthogalum Lallem. 

Zizia W. D. J. Koch 


Apium Clade” 
Ammi L. 
Anethum L. 
Apium L. 
Deverra DC. 
Foeniculum Hill 
Naufraga Constance & Cannon 
Petroselinum Hill 
Ridolfia Moris 


Arracacia Clade* 
Arracacia Bancroft 
Coaxana J. M. Coult. & Rose 
Coulterophytum B. L. Rob. 
Dahliaphyllum Constance & Breedlove 
Donnellsmithia J. M. Coult. & Rose 
Enantiophylla J. M. Coult. & Rose 
Mathiasella Constance & C. Hitch. 
Myrrhidendron J. M. Coult. & Rose 
Prionosciadium S. Watson 
Rhodosciadium S. Watson 


Conioselinum Clade? 
Conioselinum Hoffm. (in part) 
Ligusticum L. (in part) 


Heracleum Clade* 
Heracleum L. 
Malabaila Hoffm. 
Pastinaca L. 
Tetrataenium (DC.) Manden. 
Tordylium L. 
Zosima Hoffm. 


Komarovia Clade‘ 
Hansenia Turcz. 
Komarovia Korovin 
Parasilaus Leute 
Physospermopsis H. Wolff 
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TABLE 3. (Continued.) 


Pimpinella Clade® Bubon L. Registaniella Rech. f. 
Aphanopleura Boiss. Pimpinella L. 
Arafoe Pimenov & Lavrova Psammogeton Edgew. 


2 


a 


o 


EN 


ag 


The Angelica clade sensu stricto (Downie et al., 2000c). If this clade, as currently circumscribed, were to be treated at the tribal level, the 


earliest name Selineae Spreng. would apply. However, we have yet to examine material of Selinum carvifolia (L.) L., the type of the genus, 
and until we confirm its placement we refrain from formally recognizing this clade. 

The Apium clade sensu stricto (Downie et al., 2000c). This clade, if recognized at the tribal level, will be called Apieae [ined.]. 

We name this clade after the largest and earliest described genus, Arracacia. 

Comprising two species each of Conioselinum and Ligusticum but excluding the types of each of these genera (C. tataricum Hoffm. and L. 
scoticum L.). However, the type of Kreidion Raf. (K. chinensis (L.) Raf.), a synonym of Conioselinum, is included in this clade (Downie 
et al., 2000c). 

With the addition of Pastinaca, this group coincides with Drude’s Peucedaneae subtribe Tordyliinae (or tribe Tordylieae W. D. J. Koch). While 
ITS studies support strongly the monophyly of the group, some studies of plastid DNA remove Tordylium from the clade. 

Members of this clade are allied closely to tribe Pleurospermeae and the genus Erigenia ( Valiejo-Roman et al., 1998; Katz-Downie et al., 1999; 
Pimenov et al., 1999). As Pleurospermum is taxonomically complex and surrounded by many small genera of dubious affinity, we refrain from 
naming this group until these genera are examined. If Erigenia is indeed to be included in the Komarovia clade (and if the latter retains its 
present complement of genera), the tribal name Erigenieae Rydb. would apply. 

If further investigation supports this group, and if it is to be recognized at the tribal level, the earliest name Pimpinelleae Spreng. should be 
applied, as the type of Pimpinella (P. saxifraga L.) allies with this group. Registaniella is included here based on its morphological similarity 
to Psammogeton. 
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and Smyrnieae, with Scandiceae comprising subtribes Daucinae Dumort., Scandicinae 
Tausch and Torilidinae Dumort. (Downie et al., 2000a,b,c). Tribes Aciphylleae, 
Heteromorpheae and  Pleurospermeae are recently described, whereas tribes 
Bupleureae, Careae, Oenantheae, Pyramidoptereae and Smyrnieae, while maintaining 
long-standing names, are radically different in generic composition from those groups 
originally recognized. Tribe Echinophoreae is maintained as circumscribed by Drude, 
although three of its six genera have yet to be included in any molecular study, and 
the remaining genera show affinities with the Socotran Island endemics 
Nirarathamnos and Rughidia (Downie et al., 2000a,c). Whether these two genera are 
included in an expanded Echinophoreae or are treated as a separate taxon remains 
to be seen. Scandiceae subtribes Daucinae and Torilidinae approximate tribe 
Caucalideae Spreng., the subject of intensive multidisciplinary systematic investi- 
gations during the 1960s and 1970s (summarized in Heywood, 1978a). However, 
representatives of Drude's tribe Laserpitieae Benth. occur within subtribe Daucinae, 
and the genera Aphanopleura, Kozlovia and Psammogeton, treated previously in tribe 
Caucalideae (Heywood, 1978b), find affinities elsewhere ( Katz-Downie et al., 1999; 
Downie et al., 20002). 


Other major clades within Apioideae 


In addition to the aforementioned monophyletic tribes, seven other major clades 
have been recognized within subfamily Apioideae on the basis of phylogenetic analy- 
sis of molecular data (Downie et al., 1998, 2000a,b,c). These include the Angelica, 
Apium, Arracacia, Conioselinum, Heracleum, Komarovia and Pimpinella clades, and 
the generic composition of each is provided in Table 3. The Angelica and Apium 
clades were first established on the basis of cpDNA evidence (Plunkett et al., 1996b; 
Downie et al., 1998), but increased sampling has since obscured their boundaries. 
Moreover, the results of the ITS-based investigations have never been fully compat- 
ible with those of the cpDNA studies, with members of the Apium clade forming as 
many as five separate lineages basal to the Angelica clade (Downie et al., 1998; Katz- 
Downie et al., 1999). Therefore, to reconcile these groups with recent data, a more 
restricted view of these clades has been put forth (Downie et al., 2000c). As such, 
congruence is achieved in the composition of each of these clades and the bootstrap 
values supporting them are much higher. The phylogenetic position of the Arracacia 
clade is not clear, for in several studies it arises from within the Angelica clade 
(Plunkett et al., 1996b; Downie et al., 1998) and in some trees its monophyly is not 
very well supported (Downie et al., 1998). The Arracacia clade may eventually be 
subsumed within the Angelica clade. The Heracleum clade coincides with Drude's 
Peucedaneae subtribe Tordyliinae (= Tordylieae W. D. J. Koch) and its separation 
from other peucedanoid taxa was confirmed previously using seed protein immuno- 
chemistry (Shneyer et al., 1995). A strong serological similarity is also apparent 
between Komarovia and Parasilaus, two members of the Komarovia clade (Pimenov 
et al., 1999). 


TRIBES AND CLADES WITHIN APIACEAE SUBFAMILY APIOIDEAE 313 





Apioid tribes Careae, Echinophoreae and Pyramidoptereae, along with clades 
Angelica, Apium, Arracacia, Heracleum and Pimpinella (plus a number of genera of 
uncertain clade placement, see below) comprise a large assemblage of taxa termed 
the ‘apioid superclade’ (Plunkett & Downie, 1999). This superclade, confirmed as 
monophyletic in all analyses of molecular data to date and often supported by 
bootstrap values >90%, comprises the previously designated Aegopodium, Angelica, 
Apium and Crithmum clades (Plunkett et al., 1996b; Downie et al., 1998). The 
Aegopodium and Crithmum clades are now recognized as tribes Careae and 
Pyramidoptereae, respectively (Downie et al., 2000c). Previously, both Careae and 
Pyramidoptereae were considered part of an expanded Apium clade (Plunkett et al., 
1996b) or Aegopodium clade (Plunkett & Downie, 1999; Downie et al., 2000a). While 
the apioid superclade is strongly supported as monophyletic, resolution of relation- 
ships among its constituent taxa is poor. However, a specific class of chloroplast 
genome structural rearrangement has been uncovered within the superclade and 
serves to strengthen support for several lineages delimited on the basis of DNA 
sequence data. Here, at least one expansion and seven different contractions of the 
junction between the large single-copy and inverted repeat (IR) regions of the chloro- 
plast genome have been detected, each ranging in size from 1 to 16kb pairs (Plunkett 
& Downie, 2000). The frequency and large size of these IR junction shifts are unpre- 
cedented among angiosperms. Continued investigations on the phylogenetic utility 
of this molecular character in subfamily Apioideae are in order. 


Genera of uncertain phylogenetic placement 


As a consequence of the narrower circumscriptions of the Aegopodium, Angelica and 
Apium clades (Downie et al., 2000c), some genera previously placed within each are 
left without assignment to a specific group (yet they are still retained within the 
apioid superclade). These include Bifora, Conium, Coriandrum, Ferula, Levisticum 
and Prangos (and other genera identified by asterisks in Table 4). Their inclusion 
would result primarily in the blurring of the boundaries of the Angelica and Apium 
clades. The genera Cnidium, Conioselinum, Ferula, Ligusticum and Trachyspermum, 
however, are each not monophyletic and have some members assigned to specific 
tribes (Table 2) or informally recognized clades (Table 3). The Angelica clade is a 
large group that contains many genera of exceptional taxonomic difficulty (such as 
Angelica, Lomatium, Peucedanum, Selinum and Seseli). Given that many large apioid 
genera are not monophyletic (Pimenov & Leonov, 1993; Valiejo-Roman et al., 1998; 
Downie et al., 2000c), the placement of any one genus into the Angelica clade (or, 
for that matter, any other clade heretofore recognized ) must be treated as provisional 
until such genus is studied in detail. Considering those genera basal to the apioid 
superclade, Annesorhiza falls adjacent to Heteromorpheae, and Diplolophium and 
Erigenia fall adjacent to Komarovia (Downie & Katz-Downie, 1999). 
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TABLE 4. Genera of Apiaceae subfamily Apioideae included in molecular systematic studies 
but of uncertain tribal or clade placement. Asterisks (*) denote genera belonging to the apioid 
superclade. Genera designated as ‘in part’ have some members assigned to tribes (Table 2) 
or major clades (Table 3) 


Annesorhiza Cham. & Schltdl. Krubera Hoffm.* 

Azilia Hedge & Lamond* Levisticum Hill* 

Bifora Hoffm.* Ligusticum L. (in part) 
Cenolophium W. D. J. Koch* Lithosciadium Turcz.* 
Cnidium Cuss. (in part)* Meum Hill* 

Conioselinum Hoffm. (in part)* Opopanax W. D. J. Koch* 
Conium L.* Prangos Lindl.* 
Coriandrum L.* Smyrniopsis Boiss.* 
Diplolophium Turcz. Sphaenolobium Pimenov* 
Erigenia Nutt. Thysselinum Hill* 

Ferula L. (in part)* Trachyspermum Link (in part)* 


Ferulago W. D. J. Koch* 


PHYLOGENETIC UriLITY OF MOLECULAR MARKERS 


We compare the utility of seven molecular data sets used for phylogenetic inference 
in subfamily Apioideae: restriction sites, ITS, rp/16 intron, rps16 intron, rpoC1 intron, 
matK, and rbcL. While comparisons between different molecular data sets were con- 
sidered in earlier papers (e.g. Plunkett et al., 1997; Downie et al., 1998, 2000b; Plunkett 
& Downie, 1999; Lee & Downie, 2000), this study is the first to offer a comparison 
among all seven markers. The 24—33 species considered in each data set (Table 5) 
represent many of the major lineages described above and, as far as possible, were 
included in multiple studies. However, it wasn’t possible to make the data sets equi- 
table in species composition. The rbcL matrix was constructed from four separate 
studies (Olmstead et al., 1992; Chase et al., 1993; Kondo et al., 1996; Plunkett et al., 
1996a), whereas the matK study (Plunkett et al., 1996b) was carried out independently 
from those of the introns, ITS, and restriction sites. Of the 33 species included in 
each of the three intron studies, 30 species were shared. Here, Ferula, Anthriscus, and 
Aralia L. were each represented by one of two species; each species pair, however, is 
monophyletic as revealed by more comprehensive analyses (Downie et al., 2000a,b,c). 
Ten species (18 genera) were shared among the three intron and rbcL matrices, and 
10 species (17 genera) were shared among the intron and matK data sets. Sixteen 
species (19 genera) were common to both rbcL and matK analyses. The ITS study 
(24 species) excluded members of Araliaceae and Apiaceae subfamilies 
Hydrocotyloideae and Saniculoideae, as well as the basal apioids Anginon, Bupleurum 
and Heteromorpha, owing to the difficulty in aligning many of their sequences with 
those of other Apioideae (Downie & Katz-Downie, 1996; Downie et al., 2000c). The 
restriction site study (27 species) omitted Heracleum lanatum, Komarovia anisosperma, 
Physospermum cornubiense, Tetrataenium rigens, Eremocharis fruticosa Phil. and 
Petagnaea saniculifolia Guss. because data were not available (Plunkett & Downie, 
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TABLE 5. Species included in separate or combined analyses of DNA data or supertree 
construction. Treatments of Apiaceae subfamily Apioideae based on Pimenov & Leonov (1993) 
and the results of previously published molecular phylogenetic studies (Tables 2, 3). The 
symbol ‘???? denotes species of uncertain tribe or clade assignment (Table 4). Voucher and 
source information, and GenBank accession numbers for DNA sequences, are provided in 
footnote references 


Tribe (Pimenov Tribe, clade (molecular 
Species & Leonov, 1993) phylogenetic studies) 
APIACEAE SUBFAMILY APIOIDEAE 
Aciphylla aurea W. R. B. Oliv.‘ Apieae Aciphylleae 
Aegopodium podagraria L.^8 Apieae Careae 
Aethusa cynapium L.* 54 Apieae Angelica clade 
Anethum graveolens L.*5-ed4ef Peucedaneae Apium clade 
Angelica archangelica L.>®™®4° Angeliceae Angelica clade 
Angelica lucida L.^8 Angeliceae Angelica clade 
Anginon rugosum ( Thunb.) Apieae Heteromorpheae 
Raf.^e4ef.g 
Anthriscus caucalis M. Bieb.^ Scandiceae Scandiceae, Scandicinae 
Anthriscus cerefolium (L.) Scandiceae Scandiceae, Scandicinae 
Hoffm.*^** 
Anthriscus sylvestris (L.) Hoffm." Scandiceae Scandiceae, Scandicinae 
Apium graveolens L.*:5e4.ef& Apieae Apium clade 
Arracacia aegopodioides (Kunth) Smyrnieae Arracacia clade 
J. M. Coult. & Rose" 
Bupleurum chinense DC.* Apieae Bupleureae 
Bupleurum falcatum L.f8 Apieae Bupleureae 
Bupleurum fruticosum L.? Apieae Bupleureae 
Bupleurum ranunculoides L.*9 Apieae Bupleureae 
Carum carvi L.f Apieae Careae 
Caucalis platycarpos L.>®™®%° Caucalideae Scandiceae, Torilidinae 
Cicuta virosa L.*'5e4.es Apieae Oenantheae 
Conium maculatum L.*5*4& Smyrnieae ?m 
Coriandrum sativum L.^* Coriandreae 2») 
Crithmum maritimum L .* 54e Apieae Pyramidoptereae 
Cryptotaenia japonica Apieae Oenantheae 
Hassk.* edes 
Daucus carota L.*59$efs Caucalideae Scandiceae, Daucinae 
Donnellsmithia cordata Smyrnieae Arracacia clade 
(J. M. Coult. & Rose) 
Mathias & Constance? 
Enantiophylla heydeana Angeliceae Arracacia clade 
J. M. Coult. & Rose’ 
Endressia castellana Coincy*® Apieae Angelica clade 
Ferula assa-foetida L.** Peucedaneae ” 
Ferula communis L.* Peucedaneae ” 
Ferula tingitana L.54 Peucedaneae ??? 
Foeniculum vulgare Mill *-54.e.£2 Apieae Apium clade 


Heracleum dulce Fisch.® Tordylieae Heracleum clade 
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TABLE 5. (Continued.) 


Tribe (Pimenov Tribe, clade (molecular 


Species & Leonov, 1993) phylogenetic studies) 
Heracleum lanatum Michx.5d4f Tordylieae Heracleum clade 
Heteromorpha arborescens Apieae Heteromorpheae 

(Spreng.) Cham. & 
Schltd].*:^4.« 
Heteromorpha trifoliata ( Wendl.) Apieae Heteromorpheae 
Eckl. & Zeyh.^* 
Komarovia anisosperma incertae sedis Komarovia clade 
Korovin**4* 
Laserpitium latifolium L.f Laserpitieae Scandiceae, Daucinae 
Levisticum officinale Apieae ” 
W. D. J. Koch f 
Ligusticum scoticum L.*'*e4.e Apieae ” 
Oenanthe fistulosa L.* Apieae Oenantheae 
Oenanthe javanica DC.® Apieae Oenantheae 
Oenanthe pimpinelloides L4. Apieae Oenantheae 
Oenanthe sarmentosa J. Presl Apieae Oenantheae 
Osmorhiza aristata (Thunb. ex Scandiceae Scandiceae, Scandicinae 
Murr.) Makino & Y. Yabe? 
Pastinaca sativa L242! Tordylieae Heracleum clade 
Physospermum cornubiense (L.) Smyrnieae Pleurospermeae 
DC.b5e4e 
Pimpinella saxifraga L.'* Apieae Pimpinella clade 
Pleurospermum camtschaticum Smyrnieae Pleurospermeae 
Hoffm.? 
Rhodosciadium nudicaule Peucedaneae Arracacia clade 
(J. M. Coult. & Rose) Drude' 
Scandix pecten-veneris L.* 54 Scandiceae Scandiceae, Scandicinae 
Sium latifolium L.* 5e Apieae Oenantheae 
Sium serra (Franch. & Savat.) Apieae Oenantheae 
Kitag.? 
Sium suave Walt.’ Apieae Oenantheae 
Smyrnium olusatrum L.* 54 Smyrnieae Smyrnieae 
Spermolepis echinata (Nutt.) Apieae Angelica clade 
Hellerf 
Tetrataenium rigens (DC.) Tordylieae Heracleum clade 
Manden.5*4* 
Tordylium aegyptiacum (L.) Tordylieae Heracleum clade 
Lam.?5*4«e 
Torilis japonica (Houtt.) DC.? Caucalideae Scandiceae, Torilidinae 


APIACEAE SUBFAMILY 
HYDROCOTYLOIDEAE 
Bolax gummifera (Lam.) 
Spreng.*^4* 
Centella asiatica (L.) Urb.**4* 


Centella erecta (L. f.) Fern.? 
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TABLE 5. (Continued.) 


Tribe (Pimenov Tribe, clade (molecular 
Species & Leonov, 1993) phylogenetic studies) 


Eremocharis fruticosa Phi].^ 4^8 
Hydrocotyle bowlesioides 
Mathias & Constance? 
Micropleura renifolia Lag.* 
Spananthe paniculata Jacq.® 


APIACEAE SUBFAMILY SANICULOIDEAE 
Eryngium bourgattii Gouan."® 
Eryngium cervantesii 
F. Delaroche* 
Lagoecia cuminoides L.f Lagoecieae Pyramidoptereae 
Petagnaea saniculifolia Guss. *9t 
Sanicula canadensis L.°4* 
Sanicula gregari Bickn.^? 


ARALIACEAE 
Aralia chinensis L.4 
Aralia spinosa L.***:$ 
Kalopanax pictus Nakai^? 


* cpDNA restriction site study (Plunkett & Downie, 1999). 

> ITS study (Downie & Katz-Downie, 1996; Downie et al., 1998, 2000c). 

rpl16 intron study (Downie et al., 2000b). 

rpsl6 intron study (Downie et al., 2000c). 

* rpoCl intron study (Downie et al., 1998, 2000b). 

f matK study (Plunkett et al., 1996b). 

* rbcL study (Olmstead et al., 1992; Chase et al., 1993; Kondo et al., 1996; Plunkett et al., 
1996a). 


c 
d 


1999); Sanicula canadensis L. was also excluded but replaced with the closely allied 
Eryngium cervantesii F. Delaroche. The three intron, ITS, and restriction site data 
matrices shared 17 species and 20 genera (the groups of species representing Oenanthe 
and Bupleurum are also each monophyletic; Downie & Katz-Downie, 1999, and 
unpublished data), but upon consideration of rbcL and matK, the number of shared 
taxa among the seven data sets decreased substantially. Analyses of combined data 
therefore excluded data sets rbcL and matK. To reduce the number of terminals in 
the phylogenetic analysis of combined data or combined trees (discussed below), 
monophyletic groups of species belonging to the same genus (1.e. Ferula, Anthriscus, 
Aralia, Oenanthe and Bupleurum) were scored as a single taxon. 


Experimental 


Details of methodology and information on source of material, voucher deposition, and 
GenBank accession numbers are provided elsewhere (Plunkett et al., 1996a,b; Downie & Katz- 
Downie, 1996, 1999; Downie et al., 1996, 1998, 2000b,c; Kondo et al., 1996; Plunkett & 
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Downie, 1999). For the ITS study, complete ITS1 and ITS2 regions were used. Each of the 
three intron studies included data from the entire intron region, whereas only three quarters 
of the matK gene was sequenced (Plunkett et al., 1996b). Due to missing data at both 5' and 
3’ ends of rbcL for some species (Kondo et al., 1996), only 1382 of a total of 1428 positions 
were included in the multiple alignment of these sequences. For the restriction site study, data 
were obtained for 14 enzymes and only variable sites were scored (Plunkett & Downie, 1999). 
The data sets analysed herein are available upon request. 


Phylogenetic analysis 


All phylogenetic analyses were carried out using maximum parsimony, implemented using 
PAUP* version 4.0b4a (Swofford, 1998). Shortest trees for each of the eight data sets (includ- 
ing that of combined rp/16, rps16, and rpoCl intron data) were obtained using 100 random 
addition replicate searches, with tree bisection-reconnection (TBR) branch swapping and 
ACCTRAN optimization. Bootstrap values (Felsenstein, 1985) were calculated from 100 
replicate analyses, simple addition sequence of taxa, and TBR branch swapping. Pair-wise 
nucleotide differences of unambiguously aligned nucleotide positions were determined using 
the distance matrix option in PAUP* (with no provision made to account for multiple hits), 
with alignment gaps in any one sequence treated as missing for all taxa. To examine the extent 
of conflict among the three intron data sets prior to combining data, we conducted the 
incongruence length difference (ILD) test of Farris et al. (1995), as implemented by PAUP* 
using simple addition sequence and TBR branch swapping. Five hundred replicates were 
considered for each partition. The results of this partition homogeneity test revealed that the 
three intron matrices do not yield significantly different phylogenetic estimates and, as such, 
an analysis of combined data is appropriate. All trees constructed from plastid DNA data 
were rooted with Aralia (Araliaceae), with the exception of the matK trees, which were rooted 
with the closely allied araliad genus Kalopanax Miq. (Aralia was unavailable for analysis). 
ITS trees were rooted with the apioid genus Physospermum, as sequences from more distantly 
related taxa were difficult to align. 

Relationships among those genera included in the three intron, ITS, and restriction site 
data matrices were further explored upon simultaneous consideration of all data (the ‘superma- 
trix’ or ‘total evidence’ approach; Kluge & Wolf, 1993; Wiens & Reeder, 1995). Here, five 
data matrices were combined into a single larger matrix and analysed using maximum parsi- 
mony, with data for those taxa not included in any single analysis scored as missing. An 
alternative approach for simultaneous consideration of all data (as well as unique taxa) is the 
construction of supertrees, where the resultant trees from each of the separate analyses are 
combined rather than the original data used to construct them (Sanderson et al., 1998). Using 
the method of matrix representation with parsimony (MRP; Ragan, 1992), the topologies of 
each of the five source (individual or consensus) trees obtained from maximum parsimony 
analyses of the five original data sets were translated manually into data matrices of binary 
characters, with question marks used for taxa not included in a given tree (Baum, 1992; 
Sanderson et al., 1998). The binary characters from all source trees were combined into one 
matrix and analysed using maximum parsimony, as detailed above for the separate analyses. 
Bootstrap values, however, were calculated for 1000 replicate analyses; decay values were also 
obtained (Bremer, 1988). The resultant maximally parsimonious trees, or supertrees, were 
rooted with Aralia, having all of its states scored as ‘0’. We reiterate that the original data 
matrices or source trees do not share identical termini and that monophyletic groups of species 
in the same genus were treated as single terminals in these combined analyses. Therefore, 
while we present a strict consensus of supertrees, it is not a ‘strict supertree’ (Sanderson et al., 
1998) in the sense that it agrees with all the trees from which it was derived. 
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Comparative utility of molecular data 


Characteristics of the eight data sets are presented in Table 6 and the single or strict 
consensus tree resulting from maximum parsimony analysis of each data set 1s pre- 
sented in Fig. 2. In almost all instances, the phylogenetic results obtained from each 
of these reduced analyses are consistent with those of our previously published results 
when many more species were considered. The few areas of discord between the 
reduced and original analyses — most notable with regard to the rp/16 intron — are 
generally characterized by poorly supported nodes and, when collapsed, yield top- 
ologies of a similar nature. It is known that differences in species sampling may have 
an influence on the resulting phylogenetic hypotheses (Lecointre et al., 1993; Hillis, 
1996, 1998), and our results support this conclusion. 

Of the six coding (rbcL and matK) and non-coding (ITS and introns) loci 
sequenced, the ITS region is evolving most rapidly, as evidenced by the greater 
percentage of sites that are potentially parsimony informative and its higher rate of 
sequence change (that in some pair-wise comparisons exceeds 30% of nucleotides; 
Table 6). This high rate of nucleotide substitution, however, precludes the ITS 
region for deep-level analysis, as available sequence data for Hydrocotyloideae, 
Saniculoideae, and some basal Apioideae could not be aligned unambiguously with 
data from other apioid taxa. In contrast, rbcL is evolving most conservatively, with 


TABLE 6. Characteristics of eight data sets (A-H) used in maximum parsimony analyses 
of selected representatives of Apiaceae subfamily Apioideae and relatives. Only variable charac- 
ters were included in the restriction site study. The single or strict consensus tree resulting 
from each of these analyses is presented in Figs. 2A-H, respectively 








A. B. C. D. E. F. G. H. 
Restriction ITS rpll6 rpsló | rpoCl matk rbcL all 3 
sites intron intron intron introns 
No. of taxa 27 24 33 33 33 33 33 33 
No. of total characters 7718 471 1152 1083 884 1107 1382 3119 
No. of eliminated characters 0 47 167 167 97 0 0 431 
No. of informative 476 198 195 149 107 161 117 451 
characters 
No. of autapomorphic 302 104 183 158 126 192 107 467 
characters 
No. of constant characters n/a 122 607 609 554 754 1158 1770 
No. of variable characters TIS 302 378 307 233 353 224 918 
% informative characters? 61 47 20 16 14 15 8 17 
% divergence (range) 3.2-26.9 2.9-33.6 0.5-14.1 0.1-14.0 0.1-11.6 0.3-10.6 0.2-5.5 0.3-13.4 
No. of minimal length 1 2 69 8 846 224 4 8 
trees 
Length of shortest tree(s) 1594 828 664 528 369 560 393 1572 
Consistency index 0.368 0.496 0.618 0.594 0.642 0.601 0.497 0.609 
(excluding uninformative 
characters) 


Retention index 0.552 0.605 0.751 0.774 0.789 0.784 0.696 0.762 


*No. of informative characters/(no. of total characters — no. of eliminated characters). 
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only 896 of all alignment positions parsimony informative and maximum pair-wise 
sequence divergence estimates just over 596 of nucleotides. Clearly, rbcL 1s best suited 
for deep-level comparisons. Intermediate in rate of evolution is matK and the three 
chloroplast introns, with the introns rp/16 and rps16 evolving slightly faster than 
matK or the intron rpoCl. The latter two regions appear to contribute similar 
information to the phylogenetic analysis (Downie et al., 1998). The most notable 
difference, however, between coding and non-coding regions is the prevalence of 
insertions and deletions (indels) in the latter. While indels can create problems in 
alignment interpretation, unambiguously aligned gaps can provide an important 
source of phylogenetic information (Downie et al., 1998; 2000b; Downie & Katz- 
Downie, 1999). The analysis of cpDNA restriction sites provided 2—4 times more 
informative characters than any single DNA region sequenced, with levels of restric- 
tion site divergence, as assessed by mean character differences, similar to that inferred 
for the ITS region. Only in the analysis of combined intron data did the numbers 
of informative characters approach those obtained through restriction site mapping, 
although this was at the cost of sequencing over 3kb of intron DNA for each species 
compared. The continued acquisition of ITS, rps16 intron, and rp/16 intron sequences 
for molecular systematic study in Apioideae seems worthwhile, given their higher 
rates of sequence evolution over other examined loci and the ease and rapidity by 
which these sequence data can be procured. 


Phylogenetic resolution 


The results of the phylogenetic analysis of each of the seven reduced data sets and 
that of combined intron data reveal a group of species previously identified as the 
‘apioid superclade’ (AS), and with the exception of the rbcL tree (Fig. 2G), bootstrap 
support for this clade is high (90-100%). Other shared clades include those recognized 
previously as tribes Scandiceae (S; bootstrap values 77-99%), Oenantheae (O; 
63—10095), and Heteromorpheae (H; 99-100%). Data from the rp/16 intron (Fig. 2C), 
however, do not resolve a monophyletic Scandiceae. Moreover, in those cpDNA- 
derived trees in which it was included, Ligusticum (represented by its type L. scoticum) 
falls within tribe Scandiceae. In some cases, this 1s an artifact of the reduced sample 





FIG. 2. Single (tree A) and strict consensus (trees B-H ) trees derived from maximum parsi- 
mony analysis of cpDNA restriction sites (A) or DNA sequence data (B-H ). Tree H represents 
a combined analysis of the three intron data sets; the monophyletic genera Ferula, Anthriscus, 
and Aralia, each represented by one of two species in the separate analyses, were each scored 
as single exemplars. Tree summary statistics, including tree length, number of minimal length 
trees, and consistency and retention indices, are presented in Table 6. Numbers above nodes 
are bootstrap estimates for 100 replicate analyses; only values >50% are indicated. Complete 
species names are indicated in Table 5. Within subfamily Apioideae, the following major clades 
of two or more species are indicated: AS, apioid superclade; S, Scandiceae; O, Oenantheae; 
H, Heteromorpheae. 
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FIG. 3. Single minimal length 3929-step tree derived from maximum parsimony analysis of 
combined DNA sequence (rp/16 intron, rps16 intron, rpoC1 intron, ITS) and restriction site 
data (CI, excluding uninformative characters 20. 482; RI=0. 639). Numbers above nodes are 
bootstrap estimates for 1000 replicate analyses; numbers below nodes are decay values for 
trees up to 20 steps longer than those most parsimonious. The Ferula, Anthriscus, Oenanthe, 
Bupleurum and Aralia species examined in each of the separate analyses were not identical; 
however, because each group of species is monophyletic (as revealed by our earlier, more 
comprehensive analyses), data for each genus were combined and treated as a single exemplar. 
Data for genera not included in any separate analysis were scored as missing data. Complete 
species names are indicated in Table 5. The bracketed tribes and clades are those listed in 
Tables 2 and 3; the abbreviations AS, S, O and H are the same as outlined in Fig. 2. 


size as its position within the clade is not supported by parsimony analyses of a much 
larger sample of taxa (cf. restriction sites (Plunkett & Downie, 1999), rp/16 intron 
(Downie et al., 2000b), rps16 intron (Lee & Downie, 2000)). In other cases, L. scoticum 
does indeed ally with Scandiceae (Kondo et al., 1996; Downie et al., 1998); its position 
here, however, cannot be reconciled on the basis of any other kind of evidence. 

The single minimal length tree derived from maximum parsimony analysis of 
combined DNA sequence (rp/16 intron, rps16 intron, rpoC1 intron, ITS) and restric- 
tion site data is presented in Fig. 3. This tree has a length of 3929 steps, a consistency 
index (CI; excluding uninformative characters) of 0.482, and a retention index (RI) 
of 0.639. Here, all clades of two or more taxa are strongly supported as monophyletic, 
as assessed by high bootstrap and decay values. The monophyletic Scandiceae is 
sister to Ligusticum scoticum. The strict consensus supertree of 16 minimal length 
138-step trees constructed by the MRP method using five source trees (trees A-E in 
Fig. 2) and analysed using maximum parsimony is presented in Fig. 4 (CI, excluding 
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FIG. 4. Strict consensus of 16 minimal length 138-step supertrees constructed by the MRP 
(matrix representation with parsimony) method using five source trees ( Fig. 2, trees A-E) and 
analysed using maximum parsimony (CI, excluding uninformative characters 20.866; RI= 
0.958). Numbers above nodes are bootstrap estimates for 1000 replicate analyses; numbers 
below nodes are decay values for trees up to 20 steps longer than those most parsimonious. 
The Ferula, Anthriscus, Oenanthe, Bupleurum and Aralia species sampled in each of the source 
trees were not equivalent, thus to reduce the number of terminals in the supertree these 
monophyletic groups of species were each scored as a single taxon. Taxa not included in any 
single source tree were coded as missing data. Complete species names are indicated in Table 5. 
The bracketed tribes and clades are those listed in Tables 2 and 3; the abbreviations AS, S, 
O and H are the same as outlined in Fig. 2. 


uninformative characters = 0.866; RI— 0.958). This tree is highly consistent to that 
derived from the analysis of combined data (Fig. 3) and upon collapse of those 
nodes supported by bootstrap values less than 5596 the trees are almost identical, 
with the only difference between them being the position of Ligusticum. 


PHYLOGENETIC RELATIONSHIPS WITHIN APIOIDEAE 


The putative relationship among the tribes and clades so far recognized in Apiaceae 
subfamily Apioideae is illustrated in Fig. 5. This phylogeny was inferred based on 
the results of combined data and trees described above (Figs 3, 4) as well as the 
results of previously published molecular analyses (Plunkett et al., 1996b; Downie 
& Katz-Downie, 1999; Plunkett & Downie, 1999; Downie et al., 1998, 2000a,b,c). 
We emphasize that at least 23 genera included in molecular studies to date (Table 4) 
do not comfortably find a home within any of these clades, although the majority 
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FIG. 5. Summary of relationships among the tribes (Table 2) and informally recognized 
clades (Table 3) of Apiaceae subfamily Apioideae as revealed by phylogenetic analyses of 
molecular data (this study; Downie & Katz-Downie, 1999; Downie et al., 1998, 2000a,b,c). 
Not all genera examined to date have been assigned to a particular group (Table 4); many 
of these genera, however, are contained within the apioid superclade. 


of these taxa fall within the apioid superclade. Their inclusion would result primarily 
in the blurring of the boundaries of the Angelica and Apium clades. We also point 
out the provisional nature of these relationships, given that only half of the approxi- 
mately 400 genera recognized in Apioideae have been considered in molecular study. 
Many of the remaining genera, however, are ascribed to the large tribe Apieae 
(Pimenov & Leonov, 1993), and while it is impossible to know their affinity until 
they are examined, a great many will likely fall within the apioid superclade. 

Heteromorpheae are sister group to all remaining Apioideae, in line with their 
woody habit, primitive wood anatomy, southern African origin, and microfossil 
evidence (Plunkett et al., 1996b; Downie & Katz-Downie, 1999). The placement of 
Heteromorpha alongside Anethum and Apium in the ITS study of Valiejo-Roman 
et al. (1998) must be considered spurious. Successively more distally branching lin- 
eages include Bupleureae, Pleurospermeae, the Komarovia clade, and Oenantheae. 
The relationships among many of the remaining clades are equivocal. While the 
apioid superclade and tribe Scandiceae can each be unambiguously circumscribed 
on the basis of molecular data, the relationship of these two clades to Aciphylleae, 
Smyrnieae, and the Conioselinum clade is not clear. 


CONCLUSIONS 


Phylogenetic analyses of molecular data provide very little support for Drude’s 
(1898) often cited system of classification of Apioideae or for other subfamilial 
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treatments that are based largely on morphology and anatomy. Drude's tribes 
Apieae, Coriandreae, Dauceae, Laserpitieae, Peucedaneae, Scandiceae, and Smyrnieae 
are each not monophyletic. Drude's tribes Dauceae, Laserpitieae, and Scandiceae, 
however, collectively form a monophyletic group (1.e. tribe Scandiceae). Of Drude's 
eight tribes, only the small morphologically distinct Echinophoreae is retained as 
monophyletic. Of Drude's ten subtribes and of those where sampling has been extens- 
ive (therefore, excluding the monotypic subtribe Silerinae Tausch and the bitypic 
subtribe Elaeoselininae Drude, both of tribe Laserpitieae, where material has so far 
been unavailable for analysis), only Scandiceae subtribe Scandicinae and Peucedaneae 
subtribe Tordyliinae show some resemblance to groups inferred by molecular study. 
Even so, the generic compositions of these clades are not identical to those proposed 
by Drude, with Molopospermum, Grammosciadium and Rhabdosciadium now 
removed from Scandicinae and Pastinaca added to Tordyliinae. 

The results of these molecular phylogenetic analyses have helped to more accu- 
rately delimit the family Apiaceae as well as to clarify the natural subdivisions within 
its largest subfamily, Apioideae. Apiaceae comprise three major lineages, coinciding 
with subfamilies Apioideae and Saniculoideae and the Azorella clade, the latter rep- 
resenting some but not all members of Drude’s subfamily Hydrocotyloideae. Based 
on taxonomic congruence among molecular data sets, ten major clades are recognized 
within Apioideae and treated at the rank of tribe, with Scandiceae comprising three 
subtribes. Each of these tribes and subtribes is strongly supported as monophyletic 
in many separate analyses of molecular data and, occasionally, by other lines of 
evidence. In naming these clades, we followed priority of validly published generic 
names (Pimenov & Constance, 1985) as regulated by the /CBN (Greuter et al., 
2000). Seven additional clades are identified, with formal recognition pending further 
study and confirmation of monophyly. Of the 195 genera examined, 23 have yet to 
be assigned (or are assigned only in part) to a specific group. Many of these genera 
fall within the apioid superclade, thus additional phylogenetic study of this large 
clade is warranted. 

While the relationships inferred on the basis of various types of molecular markers 
are generally congruent, it has not been possible to delimit many of these newly 
circumscribed tribes and clades using obvious morphological or anatomical synapo- 
morphies (Downie et al., 2000a,b). It is hoped that future studies will add to this 
information, but if not, we would have to accept that the task of reclassifying 
Apioideae at the suprageneric level is to be accomplished on the basis of molecular 
data rather than morphology. We realize that this is a contentious issue; in the 
absence of evident and diagnostic morphological characters, several of the tribes 
recognized herein are difficult to identify or distinguish from one another. Attempts 
to construct a phylogeny for the entire subfamily using traditional (ie. non- 
molecular) characters have proven difficult, for while cladistic analyses of morpho- 
logical and anatomical data offer much insight into the phylogeny of Scandiceae 
(Spalik & Downie, 2001; Spalik et al., 2001, and unpublished; Lee et al., 2001), 
Oenantheae (Downie & Spalik, unpublished), and basal Apioideae (Roux et al., 
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1978; Van Wyk et al., 1997; Van Wyk, 2001), these data appear to be of limited 
value at higher taxonomic levels (Downie et al., 2000a; Spalik & Downie, 2001). 
Such a conclusion is not surprising, given the common dissatisfaction among system- 
atic botanists with traditional diagnostic characters used in the higher-level classifi- 
cation of Apioideae and the lack of consensus of relationships apparent among 
existing taxonomic treatments. Simply, many of the tribes and subtribes now recog- 
nized in Apioideae cannot be delimited unambiguously using morphological data 
and it may be unrealistic to think that morphological synapomorphies can be found 
for those tribes and other major clades supported by molecular studies. 

Molecular systematic research on Apioideae has been collaborative, with partici- 
pation international. These studies will continue, with emphasis on the placement of 
the 200 or so genera yet to be examined (including those heretofore refractory to 
proper phylogenetic placement) and the re-examination and refinement of the various 
tribes and clades recognized herein. Current and future revision of each of these 
major groups, incorporating both molecular and morphological data, will ultimately 
lead to the production of a modern classification for the subfamily. 
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